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We classify the more common situations in which Coulomb frustrated phase 
separation can occur and review their properties. 
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1. Introduction 

A large variety of systems with competing short and long range interactions self 
organizes in domain patterns [1]. Examples range from ferromagnetic systems [2j 
[3] to diblock copolymers[4j and, at least theoretically, neutron star matter[5]. 
Inhomogeneous states display a simple set of predominant morphologies like circular 
droplets and stripes in two-dimensional (2D) systems, and layers, cylindrical rods and 
spherical droplets in three-dimensional (3D) systems. This tendency for a common 
behavior across different systems calls for simple models which neglect the specific 
details of each system and capture the general properties. 

In this work we will discuss phase separation frustrated by the 3D Coulomb 
interaction in electronic systems but similar ideas apply also to neutron star matter 
and diblock copolymers. Tendency to phase separation in electronic systems manifests 
itself by anomalies in the electronic contribution to the free energy density /g. Two 
kind of anomalies appear often in strongly correlated systems. The first situation 
is determined by a negative compressibility region. A notorious example is the 
uniform electron gas (EG) at low density [6] but this feature appears in several other 
models [a |8l[9lll0l[IIlll2l[13lll4l[15] including neutron star matter [5]. 

The other possibility is that the inverse electronic compressibility has a point with 
a Dirac-delta-like negative divergence at some density iic- This happens when the free 
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energies of two states which are separated by a barrier, cross each other leading to a 
cusp singularity. An example is also provided by the EG. Indeed numerical simulations 
show that the Wigner crystal and the uniform phases free energy cross at some 
density [TOj. The same feature appears in several models including manganites [1 71 [TS] . 

The behavior of /e can be summarized by expanding the electronic free energy 
density fe = Fe/V around a reference density Uc as fe = a\n — nc\^ with a < and 
correspond to 7 = 2 (negative compressibility region) or 7 = 1 (cusp behavior). A 
neutral system with a < is unstable toward phase separation however in a charged 
system, this tendency is frustrated by the long range Coulomb interaction. This leads 
to the formation of inhomogeneities with a characteristic size Id determined by the 
competition between long range forces and surface energy effects. 

Regarding the experiments, the first task is to determine what can be attributed 
to Coulomb frustrated phase separation and what can not. Information on the 
mechanism of the segregation can be obtained by an observation of the morphology. 
For example in magnetoresistant manganites, inhomogeneous states with large 
insulating and/or metallic clusters with very rough interfaces have been reDorted|19). 
These resemble the domains of the random field Ising model where the surface energy 
plays a role only at very short length scales. The shape of the interfaces is determined 
by the fluctuations of the local field which stabilizes one or the other of two competing 
phases and long range forces are not needed UHl HI]- 

A quite different morphology has been observed, also in manganites, by scanning 
tunneling spectroscopy in thin films. Domains have a filamentary and droplike 
metallic and/or insulating regions in the scale of tens to thousands of nanometers 
with smooth surfaces indicating strong surface energies |22j which is consistent with 
Coulomb frustrated phase separation as predicted by Nagaev several years ago|23). 
This is reinforced by the fact that the morphology is similar to that of classical 
mesodomains[I]. Although quenched disorder can have an important role altering the 
ideal periodic domain configurations that theories would predict in the clean limit, 
the strong dissimilarity with the random field Ising model domain morphology makes 
it clear that quenched disorder is not the driving force. In this regard it is interesting 
that, even in the absence of quenched disorder, the complexity of the energy landscape 
of frustrated phase separation models can make the ordered ground state unreachable 
leading to a glassy state [24], 

It is not clear at the moment which of the two mechanism (quenched disorder or 
short and long range forces) is dominant to determine each of the anomalous properties 
of complex electronic systems and more theoretical and experimental work is needed. 
Several recent experiments, however, point to an important role of the short and 
long range force competition. In the high-temperature superconductivity context, 
scanning tunneling microscopy investigations [25j have revealed an inhomogeneous 
state of electronic superconductinglike nanopatches coexisting with poorly metallic 
regions. The domains have roundlike shapes indicating that the surface energy has an 
important role which points to a Coulomb frustrated phase separation mechanism. 

Recently a strong transport anisotropy in ultra clean Sra Ru 2O7 samples have 
been observed p6j. This striking feature is consistent with the proposal of exotic 
electronic liquid phase analogue [27] to the intermediate order states of liquid crystals 
[2H] . The idea that the domains in Sra Ru 2O7 are due to frustrated phase separation 
has already been put forward 

The phenomenology of ruthanates has a strong similarity to that observed in 
GaAs heterostructures. Several magneto-transport anisotropics arise when the Fermi 
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level lies near the middle of a highly excited Landau level [3D]. It has been proposed 
[51] that in a clean two-dimensional electron gas (2DEG) in high Landau levels, a 
uniform phase would be unstable against a charge density striped phase where the 
electron density alternates between zero and full-filling. 

Evidences for inhomogeneous states have also been reported in the 2DEG at 
zero magnetic field. Using a local probe, Hani and collaborators [32l [33] have 
shown that close to the still debated Metal-Insulator transition (MIT), mesoscopic 
inhomogeneities appear. In addition thermodynamic measuraments have shown that 
close to the MIT the compressibility departs sharply from the predictions of an 
homogeneous electron gas [34] [35]. Another interesting finding is the appearance of 
negative spikes in the electronic compressibility [32| which indicate that the transition 
from the homogeneous state to the inhomogeneous state is discontinuous as found 
theoretically for a striped Coulomb frustrated phase separated state |36j. 

In cuprates, as in manganites, the situation is more complex and different 
inhomogeneities have been reported at different length scales. At a scale of lOnm 
20 lattice constants the system appears to segregate into a pseudogap or underdoped 
phase with a large gap and a superconducting phase with a smaller gap with smooth 
interfaces in between[25]. This is consistent with Coulomb frustrated phase separation 
between an underdoped pseudogap phase and an overdoped phase [37]. This scenario 
has been reproposed on the light of recent neutron scattering experiments [35] • At 
a scale of ~ 4 lattice constants the systems shows charge ordering in the form of 
stripes[39]. By considering the Cu02 plane as a strongly correlated 2D electronic gas 
coupled to the lattice, in the presence of the long-range Coulomb interaction, it has 
been proposed[9] that an incommensurate charge density wave occurs around optimal 
doping as Coulomb frustrated phase separation with hole rich and hole poor regions 
smoothly evolving intro stripes by underdoping. 

Since domains in these system have often mesoscopic scales of several lattice 
constants, general aspect of the phenomena can be analyzed by the use of coarse 
grained models. In this work we will review the two more common coarse grained 
model corresponding to two different universality classes and their differences and 
similarities. The paper is organized as follows: In Sec. [5] we review the more relevant 
models that along the years have been proposed for 7 = 1, 2 and introduce the basic 
length scales of the problem. In Sec. [3] we discuss the main properties of phase- 
separated states for weak Coulomb interaction where a unified picture can be achieved. 
In Sec. [4] we discuss the strong interaction regime where strong differences between 
7=1,2 arise and we conclude in Sec. [5] 

2. Minimal models and typical length scales 

As discussed in Sec. [T] tendency to phase separation can be sorted in two main 
electronic free energy anomalies, corresponding to a negative compressibility density 
range or a Dirac-delta-like singularity in the compressibility. Although both situations 
can be captured expanding the free energy as fe oc a\n — nc\"' with a < and 7 = 1,2, 
a full analysis of the phase separation problem, from small to high frustration, requires 
higher order terms in the free energy expansion as follows: 

fe = a\n~n,\^ +f3\n-n,\^\ (1) 

For (3 > this provides a double-well form for the free energy with minima at 
n = Uci: SnP/2 (here Sn^ = 2[|q;|/(2/3)]^/''') and a barrier between minima of height 



Coarse grained models in Coulomb-frustrated phase separation 



4 



/o — a^/(4/3). This form of the free energy assumes two symmetric phases with the 
same compressibihty. One can consider asymmetric situations by adding cubic terms. 
Such asymmetry does not change the physics substantially and will be neglected for 
simplicity. The exception is the important case of an incompressible phase coexisting 
with a compressible phase. This limiting case constitutes a different universality class 
and will not be discussed here. A treatment can be found in Ref. [40l in 3D and in 
Ref. [36] in 2D. 

We start by introducing the "f — 2 model which is defined by the following free 



energy: 



F = d^x [aAn^ + /^An^ + c|Vn (x) |-^] (2) 



+ ^ J d^xj d^x' [n (x) - n] v{x - x') [n (x') - n] 

where n is the average charge representing the density of the rigid background, 
An — n — ric and u = |x — x'|~^ is the Coulomb interaction. Finally the gradient term 
models the surface energy of smooth interfaces and is parameterized by the stiffness 
constant c. Here the charge density plays the role of a scalar order parameter in an 
analogous way as the liquid-gas transition of classical fluids. Pattern formation within 
this model has been considered by Schmalian and Wolvnes[24j in three-dimensional 
systems, and more recently by Muratov [41] who considered the case of D-dimensional 
systems subject to a D-dimensional Coulomb interaction in the sharp-interface limit. 
This model is also relevant to describe segregation in block copolymers^. 

Since Eq. ([3]) has several parameters, to study the phase diagram it is convenient 
to measure lengths in units of ^ = y^2c/a, and define a dimensionless density 
<p{x) = 2An{x) / 6n'^ and a dimensionless free energy $ = -F/(/oC^) that, apart from 
an irrelevant constant reads: 



$ =y d^a;[0(x)2-l]%|V0(x)|2 (3) 

+ 9kJdDxJ d^x' [0(x)-^«(x-x') [0(x')-0] 
with (f) — 2iji~ nc)/5n^ and Qfj a rescaled dimensionless Coulomb coupling given by: 



a 



We see that the parameter space can be reduced to only two parameters, the 
dimensionless global density (j) and the renormalized coupling Q^. In three- 
dimensional systems Q^j coincides with the dimensionless coupling introduced in 
Ref. [H]. 

The 7—1 case has been investigated with approximate treatments in three- 
im 130] and two-dimensional systems [3S] and within an exactly solvable model in the 
limit /3 = [43] and /3 ^ 0[44j. In this case frustrated phase separation is more easily 
described by adding an auxiliary field s linearly coupled to the charge and analogous 
to a Hubbard-Stratonovich variable. Two versions of the model are possible which 
lead to similar results: s can be taken as a soft or a conventional Ising spin with 
s = ±1 where the sign distinguishes the two competing phases. We take the latter 
model which is more straightforward to analyze. It consists of a ferromagnetic Ising 
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model linearly coupled to the local charge: 

= - J ^ {s,s, - 1) - \a\ ^ ,s, {N, - TV,) (4) 



+ ^ E - + ^ E (^^ - ^) "(^^ - ^j-) (^^- - ^) 

where Sj = ±1, the index i runs over the sites of a hypercubic lattice of dimension 
-D = 2,3 with lattice constant a, the iV's are dimensionless numbers of particles per 
site and A'' their average value. The soft version replaces the Ising part with a double 
well potential 

We have written the model in the lattice for clarity but we are interested on the 
continuum limit of this model with n(x) = Ni/a^ . Uniform phases correspond to a 
ferromagnetic state in s. Inserting the two possible values of s in Eq. Q one obtains 
that uniform phases are described by Eq. ([T]) with 7 = 1, i.e. two intersecting parabolas 
with minima at ±5no/2 and a crossing point at ( the full lines in Fig.[T]). In the hard 
spin case domain walls of the Ising order parameter are sharp by construction with a 
surface tension a = 2J/a^~^ thus the Ising term can be written as ctS with S the total 
surface of interface among the two phases. It is convenient to define the analogous of 
^ for the present model ^ = AafS/a^. This represents the size that inhomogeneities 
should have for the total interface energy be of the same order as the phase separation 
energy density gain q;^/(4/3). As before we measure the energy in units of ^^Q!^/(4/3), 
lengths in units of ^ and surfaces in units of In these units cr = 1 and apart 

from an irrelevant constant one obtains the following free energy functional: 

$ = S + y" d^a;[0(x) - s(x)]2 (5) 
+ 91 fd^x f d^x' [0 (x) - 0] ^x - x') [0 (x') - ^ 



where 



fD- 



As for the model Eq. Q the parameter space is determined by the two dimensionless 
parameters <j) and Q^^. Frustration can be also measured by the parameter Q^j(^ 
which, as shown below, has the meaning of ratio between the energy cost introduced 
by frustration and phase separation energy gain. This was the parameter used in 
Refs. [421 HOl [361144]. 

In the absence of the Coulomb interaction, both models are subject to ordinary 
phase separation in a range of global densities \cj)\ < 1 as determined by the Maxwell 
(or "tangent") construction shown by the short-long dashed line in Fig. [1] The 
phase-separated state is made up of macroscopic domains with constant local densities 
(j) = ±1. For Qlj. ^ 0, macroscopic phase separation is precluded since the Coulomb 
cost grows faster than the volume in the thermodynamic limit and mesoscale domains 
appear. 

The typical size of the domains Id can be obtained by dimensional analysis. We 
define Id = ld/£.- Taking <j> ^ 1 the Coulomb cost per domain can be estimated making 
the integrals in a volume of order if as Q'jild^~^- The surface energy goes as if^^- 
Both quantities are optimized when the inhomogeneities have the size defined by 

1? = UQl- (6) 
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Figure 1. Behavior of tlie free energy densities for the electronic uniform phases 
for the 7 = 1 model Eq. (jsjl and D = 2. The dashed lines are the exact energy 
of stripe inhomogeneities for different values of the frustrating parameter. The 
short-long dashed line corresponds to the Maxwell construction (Q|j = 0). 



Another important length scale is the screening length of the Coulomb interaction 
that can be defined for two- and three-dimensional systems as: 

l^-'=[2^-\Q'n]-' (7) 

where k is a characteristic electronic compressibility of the competing homogeneous 
phases: 

K^{2pr' (7-1) (8) 
«=(2|a|)-i (7 = 2) (9) 

For both the models presented above the characteristic screening length in units of f 
can be defined as 

ff-i ^ 1/Ql, (10) 

At weak frustration Qj^ << 1 we have the following hierarchy of scales[lT] [c.f. 
Eq. ®,(lini)]: 

ls»ld»^- (11) 

This separation of lengths allows for a unified treatment of the frustrated phase 
separation mechanism at weak frustration that will be discussed in the following 
section. 



3. The weak frustration regime 

In the weak frustration regime, the effect of long-range forces can be considered as 
a small perturbation upon the ordinary phase separation mechanism. Thus, mixed 
states are expected to appear with local densities close to the two minima of the 
double-well. For systems with ^ — 2, frustrated phase separation can be analyzed 
by expanding quadratically the free energy around the two densities (f> = ±1. Then 
the bulk free energy becomes the same as for the 7 = 1 model. In addition, the 
hierarchy of length scales Eq. pTjl for Qj^ << 1 allows to consider the interface as 
sharp. A surface tension can be defined by computing the excess energy of an isolated 
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Figure 2. The left (right) panel shows the function u{u) in two- (three-) 
dimensional systems for stripes and disks (layers, rods and droplets). For the 
disk, droplet and rod inhomogeneities, the mixing energy has been computed 
referring to B-phase inhomogeneities for ^ < \ while A-phase droplets and rods 
has been considered for u > ^. The inset in the right panel shows the difference 
in u with respect to the droplet geometry to resolve the crossings between the 
different morphologies. 



interface fir. At this point, the two models become equivalent. In the rest we present 
the analysis of the 7 = 1 model which can be solved analytically for stripes in 13 = 2 
and layers in D — 3. 

Before discussing the exact solution, it is convenient to discuss an approximate 
solution which gives essentially the same result as the exact solution in weak coupling 
but in addition allows to consider other geometries. Since the state is expected 
to be similar to the macroscopically phase separated state one can take a uniform 
density approximation (UDA) in which the local density inside the domains is assumed 
constant [421 lEl HOI [Ml Si]- Comparison with the exact result shows that this gives 
very accurate results in two- and three-dimensional systems [44] . 

The low (high) density phase will be termed A {B). Defining / = ///o the free 
energy density can be put as: 

/ = (1 - V) fA + vfB {(t>B) + (12) 

where Ja/b = (0 i 1)^ and v indicates the volume fraction of the i?-phasc domains. 
&mix represents the additional energetic cost to form inhomogeneities due to the long- 
range Coulomb interaction and the interface boundary energy. For a given geometry it 
is determined by adding the Coulomb cost and the surface energy cost and optimizing 
with respect to the dimension of the inhomogeneities. It can be cast as: 

e-mix = Qfl (0B - 0a)^^"° u{v) (13) 

The dependence of the mixed states (MS) free energy upon the morphology of the 
domains is stored in the function u(t^). In Fig. [5] we plot uiy) for the different 
geometries. These functions have been evaluated in Refs. [111[35| and appear also in the 
theory of diblock copolymers [4]. The advantage of the UDA approximation is that the 
mixing energy, that encloses long-range force effects, can be computed independently 
from the specific modeling of the homogeneous phases free energy. These functions 
are valid whatever form one choses for Ja/b- 

As shown in Fig. [2] droplet-like or disk inhomogeneities are preferred for low 
volume fractions of the minority phase [v ^ Q ov v ^ 1). On the contrary stripes and 
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layers appear for ~ i. At intermediate values of v in three dimensional systems 
one finds cylindrical rods. In the diblock copolymer context the control parameter 
is the volume fraction i>. Here instead the control parameter is the global density cf). 
However for << 1, the volume fraction increases nearly linearly with the global 
density and the two parameters are practically equivalent. From Fig. [2] one sees that 
there will be a series of morphological transitions that connect droplet states near 
the homogeneous-MS transitions to the striped mixed state at <?!) ~ 0. In the present 
approximation they appear as first order however consideration of more complicated 
geometries and charge relaxation effects can change this to a smooth evolution of 
inhomogeneities that could also include "fingering" and elongation of the domains as 
in classical systems [1]. The stripe and rod phases are expected to behave as a glass in 
quench experiments [24j with labyrinth-like patterns. 

Since we are measuring energies in terms of the barrier height, which represents 
the characteristic energy gain from phase separation Eq. allows to give another 
physical interpretation to the frustration parameter: 

^2 / Char, mixing energy cost 
^ yChar. phase separation energy gain 

In Fig.[T]we show with dashed lines the typical behavior of mixed states (MS) free 
energies for the 7 = 1 model at different values of Qj^. These results are exact but the 
results within the UDA approximation at weak frustration are practically identical [44j. 
An interesting property of mixed states is that they present the "wrong" curvature; 
that is, the electronic compressibility d^fe/dcj)^ is negative. Generally, this does 
not imply a thermodynamic instability since the usual stability condition of positive 
compressibility must be formulated for the global neutral system thus including 
the background compressibility. Since in frustrated phase separation analysis, the 
inverse background compressibility is assumed to be an infinite positive number (the 
background density is fixed to the uniform average value (j>), it follows from this point 
of view that the system is in a stable MS. 

An important difference with ordinary phase separation resides in the behavior 
of the local densities of the domains. In unfrustrated phase separation the two phases 
have a constant density independently of the global density. In frustrated phase 
separation the local density of the domains decreases with an increase of the global 
density [421 [40]. Assuming that the Curie temperature is an increasing function of 
the local density, rather than the global density (controlled by doping), this could 
explain [12^ the puzzling maximum of the Curie temperature in the three-dimensional 
pcrovskite manganite Lai^jjCa^rMnOa at a; = 0.35 Ca doping [iS] not predicted by 
the conventional double-exchange mechanism [46l |47l |48] . 

4. The strong frustration regime 

The results presented above, fully determine the behavior of MS in the weak frustration 
regime. By increasing the renormalized Coulomb frustration in the models Eqs. ([I]),©, 
MS with local densities close to the reference density appear. In this case the behavior 
of the two models is radically different leading to two "universality" classes. 

We start analyzing the 7=1 universality class. As a first approximation we can 
invoke the UDA keeping in mind that in this case, since we are far from the Maxwell 
construction limit, the results should be taken with a pinch of salt. Indeed we will see 
that special care is needed in D = 2. 
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Figure 3. The full line in the left (right) panel represents the exact phase 
diagram of the 7 = 1 model for a smectic striped (layered) state in a two- (three- 
) dimensional system. The dashed line in the left panel represents the UDA 
approximation. A discrepancy between the two solutions is found for ~ ^ 
(inset of the left panel) since the UDA predicts a critical value of Q|j while the 
exact model shows a logarithmic singularity. In three-dimensional systems the 
UDA (not shown) gives qualitatively the same result as the exact solution. For 
**?/{ > Qfl c there is a common boundary (thin line) between the two homogeneous 
phases. 



For simplicity we fix the global density at the reference density ^ = and look 
for the energetic stability of a striped {D = 2) or layered {D = 3) mixed state. By 
symmetry v — 1/2 and the free energy density behaves as 

5f ^-m + \5^^ + qT<^I'^W' (7-1) (14) 

where 54> = (pB — 't'A- The first term represents the phase separation energy gain, 
the second term is an energetic cost due to compressibility effects and the last term 
is the UDA mixing energy [c.f. Eq. (fO|) ]. The condition of stability of mixed states 
reads 5f <Q. In D = 3 the last term is dominant at small 54) and combined with the 
linear term produces an energetic barrier between the homogeneous state {54) = 0) 
and the inhomogneneous state (54> 7^ 0). The quadratic term ensures stability for 
large 54- Clearly the transition will be first order with a critical frustration given 
by Qfl. c = 27/[4u(l/2)]'^ ^ 0.47. Apart form small numerical corrections this result 
coincides with the exact solution (c.f. Fig.[3|). The different power dependence of the 
phase separation energy gain and the mixing energy cost makes the UDA reliable. 

Another important prediction of the UDA is that domains have never all linear 
dimensions much larger than the screening length. This important results is quite 
general and can be derived by simple arguments (see below). A small difference 
between the UDA and the exact results arises exponentially close to (Qlj, ^, 4c — 0) 
where domains with Id » Is are possible. This region however is physically irrelevant 
since requires an unphysical tuning of the density. 

For D — 2 we are in a marginal situation. There is a delicate balance between 
the first and the third terms in Eq. (I14p which are of the same order. The transition 
looks second order with ^ = 1/^^(1/2)^ ~ 1-45. This result, however, is incorrect. 
In this marginal case charge relaxation introduces a correction which unbalances the 
two terms. Fortunately the model is analytically solvable [151 HI] for stripe and layered 
structures which are the expected morphologies close to </) = 0. In Fig. [3] we show the 
exact phase diagrams for these particular morphologies in two- and three-dimensional 
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Figure 4. Behavior of tlie cliarge density modulation at the reference density 
<j) = for a cut perpendicular to stripes in two-dimensional systems and layers in 
three-dimensional systems and domain width W = 20ls . 

systems. In D — 2 the transition lines diverge logarithmically at = and there is no 
upper limit for the mixed state. Therefore a direct first-order phase transition between 
homogeneous phases is not possible and there is always an intermediate phase. 

The same result has been first obtained in this region of the phase diagram by 
neglecting the compressibility term in the free energy cx (f)'^ |43 . Indeed, in the strong 
frustrating regime the typical size of 2D domains is exponentially larger than the 
screening length and thus the physics is determined by the slow power-law relaxation 
of the charge density (Fig|4|). The screening length plays the role of a short-length 
cutoff that removes the unphysical divergence of the charge density at the domain 
boundaries arising at /3 = 0. 

For a practical point of view strongly frustrated 2D mixed state can be quite 
difficult to observe. It appears in the exponentially narrow range of densities 




and therefore it requires an enormous accurate control of the density. It may be 
possible that this is achieved to some extent in ruthanates, where the control parameter 
is the magnetic field which allows for considerable fine tuning and the system is 
ultraclean. 

Another difficulty is that the inverse electronic compressibility 

is exponentially large and negative in all the mixed state stability range and negatively 
diverges with 1/2 critical exponent at the transition. As stated in Sec.[3l in frustrated 
phase separation models the background is assumed incompressible (Kb = 0) but real 
systems will have a small finite background compressibility > 0. This will lead to 
a volume instability analogous to the volume instability of Cerium metal[42l I49j. 

In some cases, when there is a large separation of energy scales, the background 
may behave as nearly incompressible. This can occur in ruthanates where the relevant 
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electronic phenomena occur at temperatures below a tenth of a kelvin to be compared 
with the melting temperature of the material of the order of hundreds of kelvin. 
In the case of a 2D electron gas with a ionic background there is also a similar 
separation of energy scales. Even more, from an elastic point of view, the background 
is three dimensional and therefore behaves as practically incompressible. In field- 
effect-transistors where the background is made of mobile charges, the situation is 
quite different and the background compressibility can not be neglected |50j. 

Differences arising between two- and three-dimensional systems are mainly 
due to the different behavior of the charge density profile inside the domains. 
Phase separation energy gain stems from the region where the electronic density 
is significantly different from its average value [42]. In two-dimensional systems the 
power-law behavior of the charge density allows to gain phase separation energy even 
far from the boundaries. On the contrary, in three-dimensional systems the charge 
density decays exponentially from the domain boundary on the scale of the Thomas- 
Fermi screening length. Thus the system gains phase separation energy in a range of 
the order of the screening length from the domain boundaries. Regions far from the 
boundaries produce an exponentially small energy gain. Thus if phase separation is 
favorable the system adjust itself in such a way to eliminate all these regions. Domains 
satisfy a maximum size rule such that every point of the domain is at a distance of 
the order of the screening length or smaller from the boundaries. This rule proves to 
be quite general and is independent of the model. The maximum size rule allows for 
arbitrary large inhomogeneities in two-dimensional systems since one of the dimensions 
is already smaller than l^^ . 

Now we discuss the case of 7 = 2 at strong frustration. The separation of length 
scales discussed in Sec.|3]is no longer valid, the UDA becomes unreliable and the effect 
of smooth interfaces plays a prominent role. 

Computing the static response to an external field in momentum space one obtains 
the following charge susceptibility: 



We remind that the system is considered to be _D-dimensional but embedded in the 
usual three dimensional Coulomb interaction. The Fourier transform reads: 



X diverges for Z? = 2, 3 at a characteristic finite wavevector on an instability line 
QK^). Above g|_^ = l/(27r) in 3D and Qj^ ^ = (4V2)/(3\/37r) in 2D, the system is 
always homogeneous. On entering in the unstable region a sinusoidal charge density 
wave (SCDW) occurs. Thus inhomogeneities are quite different form the mesodomains 
of the weak coupling regime (Sec. [3]). The crossover between this two regimes will be 
discussed elsewhere [5T] . 

A related mechanism has been proposed in cuprates predicting charge ordering 
instabilities and other anomalous properties in accord with experiment [9] . In classical 
systems this is offer referred as the "microphase" separation transition [H [41] . 




1 



Xiq)= <?2-2 + 60' + 
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5. Conclusions 

In this work we reviewed the main aspects of frustrated phase separation in charged 
systems considering two kind of compressibility anomalies that generally occur in 
strongly correlated systems. For the different coarse-grained models that have been 
introduced in the literature, the outcome of long-range forces can be measured by a 
dimensionless parameter that defines the amount of frustrating effects. Frustration 
tends to reduce the range of density where a mixed state appears. Thus uniform phases 
are possible at densities where in the absence of long-range forces, phase separation 
would occur. This situation is in accord with thermodynamic measurements [Ml 
of the uniform 2DEG where one finds stability of the uniform phase with negative 
compressibility. 

In the presence of frustration the mixed state consists of domains of mesoscopic 
size with various geometries depending on the control parameters. A sequence of 
morphological transitions occurs that has a strong similarity to those in other systems 

mil- 

When frustrating effects are a small perturbation, a unified description of MS 
is allowed and a simple approximation in which the density inside the domains is 
assumed constant, gives very accurate results. 

At strong frustration, the properties of the mixed state strongly depends upon the 
particular anomaly in the compressibility and one can define two different universality 
classes. For systems with a cusp singularity in the electronic compressibility (7 — 1) 
the system dimensionality plays a key role in the charge segregation mechanism. A 
critical frustration Qr^c, above which mixed states are not possible, exist only in 
£» = 3. 

In systems with a negative electronic compressibility region (7 = 2) a critical 
value of the frustration exist Qr^c for both = 2, 3. Close to Qr^c inhomogeneities 
are SCDW. 

The maximum size rule says that domains can not have all linear dimensions much 
larger than the screening length independently of Qr. Thus a necessary condition for 
the applicability of a coarse grained treatment is that the three dimensional screening 
length must be much greater than the typical microscopic lengths like the interparticle 
distance and the lattice constant. Therefore mesoscopic domains can be expected in 
systems with small compressibility as bad metals or systems close to metal-insulator 
transitions or systems with very anisotropic electronic properties (i.e. are nearly 
insulating in one direction). Interestingly most of the materials where mesoscopic 
domains have been found have these characteristics. 
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